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1 Introduction 

The Enclosure Method introduced in [11, 10, 13] has become a well-known guiding principle in 
attacking various inverse obstacle problems [33] governed by partial differential equations. It is 
simpler than the Probe Method which has been introduced in [8, 9, 12]. 

The Enclosure Method aims at obtaining information about the geometry of unknown dis¬ 
continuity. The method consists of three steps listed below: 

• choosing a special solution v depending on a large parameter r > 0 and independent of 
the unknown discontinuity; 

• constructing a so-called indicator function of independent variable r by using observation 
data and v\ 

• studying asymptotic behaviour of the indicator function as r —> oo. 

From the asymptotic behaviour of the indication function we find a domain that encloses 
unknown discontinuity. The Enclosure Method is quite flexible and its realization depends on 

‘This paper is the author version of: Ikehata, M., New development of the enclosure method for inverse 
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L. Munteanu, V. Chiroiu), Vol. 2, Editura Academiei, Bucharest, Romania. 
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the choice of v in the Hrst step and whether the observation data in the second step depend on 
v or not. 

Now we have many applications of this flexible method to various inverse obstacle problems 
governed by elliptic partial differential equations or systems. See [ 22 ] for the systematic expla¬ 
nation of the Enclosure Method from the beginning and also [15, 19, 26, 31, 32, 35, 36] and [37] 
with references therein for further applications. For a nonlinear partial differential equation we 
cite also a recent remarkable work [3]. 

It was the paper [14] which opend the door to various possibilities of the Enclosure Method in 
the time domain inverse obstacle problems governed by the heat or wave equations in one-space 
dimension. Now we have several papers [28, 16, 29, 18, 27, 30] in which the range of application 
of the Enclosure Method has been extended to inverse obstacle problems governed by parabolic 
or hyperbolic equations over a finite time interval in three-space dimensions. 

The aim of this chapter is to make a review of the recent results using the Enclosure Method 
on inverse obstacle problems governed by the wave equation and the Maxwell system in time 
domain. We also describe some of unsolved problems related to further possibility of the En¬ 
closure Method itself. Those are not mentioned in the expository paper [22] and survey paper 

[17]- 

2 The enclosure method for inverse obstacle scattering in time 
domain 

The description of the problem is simple. Send a wave and observe the reflected wave by an 
unknown obstacle (discontinuity). What information about the obstacle can one extract from 
the observed wave? This type of problems have their origin in sonar, radar, nondestructive 
testing, etc.. 

Recently, using the Enclosure Method as the guiding principle, we considered the problem 
under the constraints: sending at most finitely many waves at a finite distance from the obstacle; 
observing a reflected wave over a finite time interval and thus at a finite distance from the 
obstacle; sending and observing places are same. In particular, we use neither the asymptotic 
behaviour of the wave as time goes to infinity nor the far field profile. In this section we present 
some of recent results from [20, 24] and their applications. 

Let D be a non-empty bounded open subset of R 3 with C' 2 -boundary such that R 3 \ D is 
connected. Let 0 < T < oo. Let / G L 2 (R 3 ) satisfy supp / fl D = 0. 

We denote by Uf the (weak) solution of the following initial boundary value problem for the 
wave equation: 

dfu — A u = 0 
c)v 

— - r y(x)dtu — (3{x)u = 0 on dDx]0, T[, 

< OU ( 1 ) 

u(x, 0) = 0 in R 3 \ D, 

dtu(x, 0) = f(x ) in R 3 \ D. 

where v denotes the unit outward normal to D on dD, (3 € L°°(dD), 7 G L°°(dD ) and j(x) > 0 
a.e. x G dD. We omit the description about the solution class taken from [5]. See [20] for the 
description. 

The role of 7 > 0 can be seen from the formal computation 

£'{t) = — [ 'y(x)\dtu\ 2 dx < 0, 

JdD 


in (R 3 \D)x ]0, T[, 
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where 


£(t) = \ j _{\dtu\ 2 + \X7u\ 2 )dx + \ l /3(x)\u\ 2 dS, t € [0, T}. 

2 Jr 3 \d 2 Jqd 

We think that the distribution of the values of 7 and (3 on 8D is a mathematical model of the 
state of the surface of the obstacle. 

Let B be the open ball centred at p with very small radius 7 and satisfy B fi D = 0. We 
denote by xb the charactersitic function of B. 

Problem 2.1. Generate u = Uf by the initial data / = xb and observe u on B over time 
interval ]0, T[. Extract information about the geometry of D , 7 and /3 from the observed data. 

The correspondence (D, 7 ,/?) '—> ^|sx]o,T[ is nonlinear. Therefore, Problem 2.1 becomes a 
nonlinear problem. 


2.1 Indicator function 


The Enclosure Method in time domain also introduces an indicator function like the Enclosure 
Method in frequency domain. It starts with introducing a special solution with a large parameter. 

In what follows we always choose / = xb unless otherwise specified and Uf is the solution 
of ( 1 ). 

Let r > 0 and v = Vf( ■, r) € H 1 ( R 3 ) be the solution of 

(A — t 2 )v + / = 0 in R 3 . (2) 


Vf has the expreesion 





p—r\x—y\ 

7- ydy. 

\x - y | 


(3) 


Define 


Wf(x,r) = / e Tt Uf(x,t)dt, x G \ D. 

Jo 


Using Vf and wj, we define the indicator function of r: 

Ib{t) = / (wf - v f )dx. 

J B 


(4) 


This indicator function looks different from the one in the previous version of the Enclosure 
Method [10]. So someone may have a question: why should it be called the indicator function? 
However, in [20], it is shown that the asymptotic formula 



is valid for an arbitrary fixed T < 00 and bounded domain Q with a smooth boundary, such 
that Dell, BnO = 0 and R 3 \ Q is connected. 

And from (4) we have the space-time expression 


Idn 


d Vf 

~du Wf 


dwf 

~dv Vf 


dS = 


IM 


d(e rt Vf) 
dv 


u f 


du 


4 (e Tt Vf ) ] dSdt, 


du 


where M = d£lx ]0, T[. Note also that: (2) implies that e~ rt Vf satisfies the wave equation in a 
neighbourhood of D. 

Theorefore, one can say that Ib{j) is essentially similar to the indicator function in the 
previous version of the Enlosure Method. It is a space-time version. 
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2.2 Qualitative state of the surface, distance and direction 

The following result is the starting point of the Enclosure Method in time domain. 

Theorem 2.1([20]). Let T > 2 dist(D,B). Let C be a positive constant. 

We have: if j(x) < 1 — C a.e. x G dD, then there exists ro > 0 such that Ib(t ) > 0 for all 
t > To; if^f(x) > 1 + C a.e. x G dD, then there exists To > 0 such that Ib(t ) < 0 for all r > ro- 
Moreover, in both cases, the formula 

lim — log \Ib(t)\ = — 2dist (D, B), (5) 

T->-00 "f 


is valid. 

Define dgu(p) = inf \x — p\. We see that knowing dist (D,B) is equivalent to knowing 
x£d D 

dgD(p) since dist (D,B) = dgo{p) — rj. Thus Theorem 2.1 yields the sphere \x — p\ = dgo{p) 
whose exterior contains D and on which at least one ponit on dD exists. Thus this should be 
called an enclosing method by using the exterior of a sphere. And also, roughly speaking, we 
can know the qualitative state of the surface of the unknown obstacle, that is whether 7 >> 1 
or 7 << 1 by the signature of indicator function Ib(j) for a large r as Theorem 2.1 states. 

Finally we present a procedure for making a decision around p whether given direction u G S 2 
the point p + dgo(p)w belongs to dD or not provided dgo(p) is known. 

Fix a large T and small s G]0, dgD(p)[- Give direction oj G S 2 choose an open ball B' centred 
at p + su> such that B' is contained in the open ball centred at p with radius dgn (p) ■ 

Step 1. Generate up by the initial data / = xb' an d observe Uf on B' over time interval ]0, T[. 
Step 2. Calculate dgo{p + su) from the data obtained in Step 1 via (5) in Theorem 2.1 in the 
case when B is replaced with B'. 

We always have dgo(p + sea) > dgo(p) — s. Moreover, it holds that: 

• if dgn{p + su) = dgoip) — s, then p + dgo(p)uj is on dD; 

• if dgjg{p + su) > dgjjip) — s, then p + dg£>{p)u is not on dD. 

Therefore one can make a decision around p whether p + dgr)(p)u is on dD or not. 


2.3 A sketch of the proof of Theorem 2.1 

The proof of Theorem 2.1 consists of three parts as described below. 
Lemma 2.1. We have, as r —> 00 


IWI WVD) = 0(e- rdist ^.Bt + e -rT )t 

= 0(T(e- rdi ’C D B) + e~ rT )) 

and 

WRWL^aD) = 0(T l '\e-^ st ^ +e- T )), 

where R = Wf — Vf. 

A brief outline of the proof of Lemma 2.1 is as follows. It follows from (1) and (2) that R 
satisfies 


(A — t 2 )R = e 


-tT ] 


in R 3 \ I). 


dR 


( dv 


dn 


\dn 


Ap — cR = — ( -A— cv) + e rT G on dD, 


( 6 ) 
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where 


c = c(x, t) = 7 {x)t + /3(x), 

F = F(x, t) = dtu{x , T) + ™(x, T), 

G = G(x) = 7 (x)u(x,T). 

( 6 ) and integration by parts give 

f _{\VR\ 2 + t 2 \R\ 2 + e~ rT FR)dx 
Jr 3 \d 

+ L{ c ' R] 2 -{r,- cv ) R+e ~ rTaR } dS = 0 ‘ 


( 7 ) 


( 8 ) 


Using a trace theorem [7] and the assumption 7 ( 2 ;) > 0 a.e. x G 51?, from ( 8 ) one can easily 
deduce the conclusion. See Lemma 2.1 and (2.28) in [20]. 

Let us continue the sketch of the proof of Theorem 2.1. Since v satisfies (2), we obtain 


/ _ fRdx = I ^-RdS+ f _(VvVR + T 2 vR)dx. 

Jr 3 \D JdD ou Jr3\D 


R 3 \D 


(9) 


On the other hand, from ( 6 ) we obtain 

0 = J |ci? — — cvj | vdS + J _(V R ■ Vv + r 2 Rv)dx 


+e 


-tT 


Ir 3 \d 

[ _Fvdx+ [ GvdS) . 
/R 3 \D JdD J 


( 10 ) 


Taking the difference of (9) from (10), we obtain 


dv 

du 


— cv ) vdS 


—e 


-tT 


/r 3 \d 


Fvdx + [ GvdS ] 
JdD J 


Then, applying ( 8 ) to the first term on this right-hand side, we obtain another expression 

I B (t)= f _(|Vi ?| 2 + T 2 | i ?| 2 )( fe + [ c\R\ 2 dS + [ (-k~ — cv) vdS 

Jr 3 \D JdD JdD \OV J 


+e 


-tT 


/r 3 yd 


FRdx+ / GRdS - 
JdD 


/R 3 \D 


Fvdx — 


[ GTciS) 
JdD J 


( 11 ) 


It is easy to see that, as r — » oo IMIl 2 (r3\;d) = ^( T 2 ) an( l H v llL 2 (aR) = 0(e T( iist (D,B)^ Thus 
from this, (7) and Lemma 2.1 we see that the last term in the right-hand side on (11) has bound 
0(T~ 1 e-' rT ) as r —> oo. Therefore we have, as r —> oo 


where 


I b {t) = E(t) + J(t) + 0(t x e rT ), 

E(t)= [ _(\VR\ 2 + r 2 \R\ 2 )dx+ [ c\R\‘ 
JR 3 \D JdD 


( 12 ) 


dS 
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and 


J(t)= L - cv ) vdS ■ 


A combination of this and Lemma 2.1 yields the following estimates on the bound of the 
indicator function. 

Lemma 2.2. We have the following asymptotic estimates: 

(i) if 0 < 7 (x) a.e. x £ dD, then as r —> oo we have 


Ib{t) > J{r) + 0{r- l e- TT ) 


(13) 


(ii) if y(x) > C a.e. x £ dD for a positive constant C', then as r —> oo we have 


Ib{t)<J(t)+[ - 

JdD c 


dv 

du 


— cv 


2 

dS + 0{T- 1 e~ TT ); 


(14) 


(in) as r 


oo we have 


\I B (r)\ = o(T 2 e~ 2Tdist(D ’ B)) +T~ 1 e- rT ). 


(15) 


A brief outline of the proof of Lemma 2.2 is as follows. From (3), (12) and Lemma 2.1 we 
have (15); (13) is clear from (12) and the positivity of E(t) for r >> 1 which is a consequence 
of the trace theorem [7]. We present here a sketch of the proof of (14). See [20] for the full 
proof. Assume that y(x) > C' a.e. x £ dD for a positive constant C'. Rewrite (8) as 



^2|VR| 2 + 2r 2 


R + 


e~ TT F 


dx 


+ 



2c 




(16) 




/ 1 

IdD 2c 



— e 


-tT 


G 


2 

dS. 


Since we have 


t 2 |R| 2 < 2r 2 


,-tT i 


R + 


2 T 2 


,-2 tT 


+ 


\Ff 


2 r 2 


and 


c|ii| 2 < 2c 


R -h{(^- Cu] -e- rT G 


+ 


1 

2c 


dv \ ~ rT r\ 2 

- - cv] - e G | 

ov 


noting the trivial inequality |Vi?| 2 < 2|Vi?| 2 , from (16) we obtain 

2 rT 

/R3\D 

Writing 


F(r) < ——j— I _\F\ 2 dx+f -(^-cv) 
T z J R3\d JdD c \au ) 


- e~ rT G 


dS. 


(17) 


dv \ _ rT 2 dv 
0 ^- cv )- e G = B~ v 


— CV 


2 -2(^-cy)e-^a + e-^\G£ 
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we have 


/ “ 

IdD C 




dS 


j i 

IdD C 



2 


dS 


+ 0 (re _T ^ dist ( Z) ’ B )+ T ) 4 - e _2rT ). 


Therefore (17) yields 


£(r) < / 

Ja 


1 

dD C 


dv 

dv 


— cv 


dS + 0(re- r(dist (AB)+T) + e - 2rT ). 


( 18 ) 


Now a combination of (18) and (12) gives (14). 

By virtue of Lemma 2.2 it suffices to study the asymptotic behaviour of two Laplace type 
integrals in (i) and (ii) from below and above, respectively. For this we have the following 
estimates. 

Lemma 2.3. Let C be a positivie constant. 

(!) If 0 < 7 (x) < 1 — C a.e. x € dD, then there exist positive numbers p, C' and tq such 
that, for all t > tq 


J(t) > c' T -^e- 2Tdis ^ D ’ B \ 


(ii) If 7 (x) > 1 + C a.e. x € dD, then there exist positive numbers p, C and tq such that, 


for all T > Tq 


J{t)+ f - 
J dD C 



2 

dS < 


C ' T -U e -2rdist(D,B) 


For the proof of Lemma 2.3 we refer the reader to [20]. The proof given therein is based on 
an argument done in [29] and covers more general /. Here we describe roughly why y(x) = 1 is 
exceptional in Lemma 2.3. 

Applying the mean value theorem [4] to (3), we have 


v(x) 


T 3 


e ~T\x~p\ 

|.T — p I 


x € R 3 \ B, 


where </?(£) = £ cosh £ — sinh £. 

Define A goip) = {q £ dD \ \q — p\ = dgn^p)}. We call A qd{p) the first reflector from p to 
dD and the points in the first reflector are called the first-reflection points, going from p to dD. 
Let x € A 8d{p)- Since u x = (p — x)/\x — p |, from the expression above we obtain 

dv ( l \ 

— = tv 1 —|- r ~ TV 

ou \ t\x — p\ J 


and hence 


dv 

dv 


— cv ~ 


t( 1 — y(x))u. 


Since all the points in A qd(p) attains the minimum of the function: dD 3 x i —> \x — p\, roughly 
speaking, one may expect, as r —> oo 


J(r) ~ r / (1 — 'y)v 2 dS 

JdD 


and 


At) + / ~ 

JdD c 


dv 

dv 


— cv 


, Q f 1 (dv \ dV r Q 

dd= - w- cv —dS ~ t 

JdD c \dv J dv 


1 T 2 


v z dS. 


IdD 7 
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These suggest (i) and (ii) of Lemma 2.3. Note also that 1 — 7 ( 2 ;) < (1 — 'y{x))/'y{x) if 7 ( 2 ;) > 0. 

Now it is easy to see that from Lemmas 2.2 and 2.3 one obtains (5) and other statements of 
Theorem 2.1. In the proof we never make use of the idea of geometrical optics which is classical. 
Everything can be done in the context of the weak solution of [5] and main tool is just integration 
by parts. Note that in Theorem 2.1 7 is just essentially bounded on dD and thus may have, for 
example, a first kind of discontinuity. 

2.4 Curvatures and counting number 

Let z £ R 3 and 0 < r. In what follows we denote by B r (z) the open ball centred at 2 and with 
radius r. 

Let q £ dD. Given v £ T q (dD) define S q (dD)v = — — (i' q ^)\t=o, where q(t) £ dD, q( 0) = q 

and dq/dt( 0) = v. We have S q {dD)v £ T q (dD). The operator S q (dD) : T q (dD) —> T q (dD) is 
called the shape operator (or Weingarten map) of dD at q derived from u. The shape operator is 
symmetric with respect to the induced inner product on T q (dD) and its eigenvalues k\(q) < k 2 (q) 
are called the principle curvatures at q. Kgjg(q) = k\(q)k 2 (q) and Hgjg(q) = (ki(q) + ^ 2 (<?))/2 
are called the Gauss and mean curvatures at q, respectively. 

Let p £ R 3 \ D. Let q 1 £ dB dgD ^(p) and S q i(dB dgD ^ (p) ) denote the shape operator of 
dBd dD (j>)(p ) at q' derived from the unit inward normal to dB dgD ^fp). If q £ A go ip), then we 
have q £ dB dgD{p) (p ), T q (dD) = T q (dB dgD ( p) (p )) and S q (dB dgD{p) (p )) - S q (dD) > 0. Since 
S q( dB d dD {p)(p)) = (1 /dg D (p))I, we have 

det (S q (dB dgD{p) (p)) - S q (dD)) = (A - h(q))( A - k 2 (q)), (19) 

where A = 1 /dgo(p). 

Theorem 2.2([24])). Let 7 = 0 . Assume that dD is C 3 and (3 £ C 2 (dD); Ago(p) is finite and 
satisfies 

det (S q (dB dgD (p ) (p) ) - S q (dD )) >0, \/q £ Ag D (p). 

If T > 2 dist(D,B), then we have 

where 

A dD(p ) = 7 

qeAenlp) yJdet(S q (dB dgD(p) (p )) - S q (dD )) 

Using Theorem 2.2, one can give a procedure for extracting the curvatures at a known first 
reflection point. More precisely, let p £ R 3 \ D and q £ Agjg(p). From (19) we have 

det (S q (dB dgD{p) (p)) - S q (dD)) = Q( A) = A 2 - 2H dD (q)X + K dD (q), (22) 

where A = 1/dgoip)- Replace p —> pj = p — SjV q , j = 1,2, 0 < si < s 2 < dgoip)- 
Then A g D (pj) = { 9 } and det (SfidB^p.^pj)) - S q (dD)) > 0 since SfidB^^pj)) > 
Sq(dB dgD ( p )(p )) and S q (dB dgL){p) (p )) - S q (dD) > 0 (q attains mm xedD \x-p\). 

Let B\ and B 2 denote two open balls cetred at p — sjv q , j = 1,2, respectively with 0 < si < 
s 2 < dgo(p) and satisfy B\ U B 2 C R 3 \ D. Let T > 2 max dist (D, Bj) and / = XBj- Applying 
( 21 ) to this case, we obtain 
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( 20 ) 

( 21 ) 



where A j = 1 /dgn(pj )■ Since dist (D,Bj) = dgo(pj) — Sj and dgo{pj) = dgr>(p) — Sj, one can 
know Q(Xj) with j = 1, 2 from Uf(x, t ) given at all (x, t ) € BjX ]0, T[ for / = XBj with j = 1, 2. 
Then, solving the system 


/ -2Ai 1 \ 

1 Haniq) \ / 

V -2A 2 1 ) 

V K 9D (q ) J \ 


Q(Ai)-A? 


(23) 


we obtain both Kgjg(q ) and Hgr>(q). Thus, one can know an approximate shape of the obstacle 
in a neighbourhood of q. Note that if dg£>{p) —> oo, then A j —> 0 and thus it will be difficult 
to extract Hgo(q) from (23). 

Another simple corollay is a formula for counting the number of unknown spherical obstacles 
with the same and known radius nearest to the center of the support of /. Assume that D = 
B e (x i) U • • • U B e (x m ), where B e (xj), j = 1, • • •, m is the open ball centred at Xj with a known 
radius e > 0 and B e (xi) D B e (xj) = 0 if i ^ j. 

Given p G R 3 \ D it is easy to see that: A qd{p) consists of finite points; (20) is satisfied; 
there exists at most one first reffection point going from p on each dB e (xj). Therefore, one can 
apply Theorem 2.2 to this case and obtain the formula which enables us to know the counting 
number of the balls which are closest to the centre of B, that is, 


tt ^3d(p) = 


\ddD(,P ) 



2 

TT 


( diam B dgD ( p ){p) 


diam B 


lim r 4 e 2Tdist ( D ’ B )/ B (r), 

T->-00 


where 


X.8d(p) = \xi + e-p —— 
l | P~Xi 


I P - Xi 


min |p 

j 



2.5 A sketch of the proof of Theorem 2.2 

Let 7 = 0. Integration by parts yields 

J{ T ) = [ (|Vu | 2 + T J \v\ 2 )dx — [ (3\v\ 2 dS. 

Jd JdD 

Applying a trace theorem [7] to the second integral on this right-hand side, we see tha J(r) > 0 
for all r >> 1. Then, we have, as r —>• oo 

E( T ) = J(t)(1 + 0(t- 1 / 2 )) (24) 

and thus from ( 12 ) we obtain 

Ib(t) = 2 <7 (t)( 1 + 0(r~ 1/2 )) + 0(r - 1 e~ rT ). (25) 

Using the Laplace method [2], one can expand J(r) under the condition (20) and we find its 
leading term which contains information about the geometry of dD at all the first reffection 
points, going from the centre of B to dD. This yields (21). 

Thus the crucial point of the proof of Theorem 2.2 is the derivation of (24). It is a combination 
of a modification of the Lax-Phillips reffection argument in [34] and a change of a dependent 
variable near dD. Here we describe the idea of the derivation of (24) in the simplest case 
7 = /? = 0. 

Since G = 0, it follows from ( 8 ) 

E(t)= f ^-RdS — e~ rT [ _FRdx 
JdD dv Jr.3\D 
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and applying the boundary condition in ( 6 ) to J(r), we obtain 


E(r)-J M = / 8 D (^ii + |b)dS 


— e 


-tT 


FRdx. 


IR3\D 


Choose v(x),x G R 3 \ D in such a way that v has a compact support and satisfies 

v = v on dD, 

dv dv 

= — 5 - on dD - 

Of Of 

Then, integration by parts and ( 6 ) gives 


( 26 ) 


J dD 

Hence we obtain 


(dv dR. 

R~ R+ 7T v ) dS = 
\ou a v 


/ _(A — t 2 )v ■ Rdx — e rT / _vFdx. 

'r 3 \d Jr 3 \d 


E(t) — J(t)= I _(A — t 2 )v ■ Rdx — e tT f _F(R + v)dx. 
Jr 3 \d Jr 3 \d 


r 3 \D 


(27) 


The point is the choice of v. Let x r denote the reflection in the tubular neighbourhood {x G 
R 3 \ D | dgoix) < 2c>o} of dD with sufficiently small d'o > 0. It is given by x r = 2 q(x) — x, 
where q(x) denote the unique point on dD such that dQ£>{x ) = \x — q(x )|. It is known that q(x) 
is C 2 for x G R 3 \ D with dgr>(x) < 25o if dD is C 3 ([ 6 ]). Choose a cutoff function (j)$ with 
0 < 5 < 5q which satisfies 0 < 4 >s(x) < 1 ; 4 >s(x) = 1 if dg£)(x) < 5 ; 4 >s(x) = 0 if dgr>(x) > 25; 
|V&(s)| < C 8 - 1 ; \V 2 Mx)\ < C5~ 2 . 

Define 

v{x) = 05 (x)v(x r ). 

Clearly (26) is satisfied with this v. A direct computation gives 


(A — t 2 )v(x) = 4 > s(x)dgD(x ) ^ aij(x)(didjv)(x r ) + (lower order terms), 


(28) 




where dij(x) with i,j = 1,2,3 are C 1 for x G R 3 \ D with dgn{x) < 25 q and independent of r, 
4>S and v. Note that the computation is based on the formula 

2 q'(x) -1 = 1 - 2 v q( . x) ® v q ( x) - 2 dg D (x)(u q{x) Y, 

where x G R 3 \ D and dgr,(x) « 1; q'(x) denotes the Jacobian matrix of the map: x \—)• q(x). 
It is a consequence of the expression q{x) = x — dgo{x)v q ^ and the formula V (dgoix)) = v q (x)- 
The point is dgjg{x) in the first term on the right-hand side of (28). By using the change of 
variable x = y r we have 


= £ 


1,3 


’ D 


/ _(A — t 2 )v ■ Rdx 

Jr 3 \d 

h{y r )dgD(y r )a ij (y r )(d i d j v)(y)R(y r )J{y)dy + 


(29) 


where J(y) denote the Jacobian of the map: y 
integration by parts yields 


y r . Since dg D (y r ) = dg D (y) = 0 on dD, 


id 


(t > s(y r )dg D {y r )ai j (y r )(d i d j v){y)R(y r )J(y)dy 


= - f d i {(f) 5 (y r )d aD (y)a ij (y r )R(y r )J(y)}d j v(y)dy. 
Jd 
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Hereafter simply estimating this right-hand side together with other terms in (29), we obtain 



— t 2 )v ■ Rdx 


<C((5||Vi?niL 2 (D,)+^ 1 p r ||L 2 (D,))l|V.;|| L2{D) +r 2 ||i?niL 2 (D,)lbllL 2 (D)), 

where D s = {y € D \ d dD (y) < 25} and R r {y ) = R{y r ). 

Here we note that 

l|Vv|| L 2 (D) < J(r) 1/2 , \\v\\ L 2 (d) < t~ 1 J(t ) 1/2 

and 

l|ViT| \ lH d s ) < CEir) 1 / 2 , \\R r \\ LHDs) < r^Eirf/ 2 . 

Choosing 5 = t^ 1 / 2 with r >> 1, we finally obtain 



t 2 )v ■ Rdx 


< Ct~ 1/2 {E(t)J{t)) 1/2 . 


From these together with (27) and the estimate ||u 


I E(r) - J(t) | < C{T- l / 2 {E{r)J{T)f/ 2 + 


L 2 (R 3\o) < C\\v\\l*( D ), we obtain 
e~ rT E(t) 1/2 + e~ rT J(t) 1/2 ) 


(30) 


and hence 

(1 - 2Cr -1 - 2e~ TT )E(r) < (1 + 2 C + 2e~ rT ) J(r) + 4e" rT . 

Therefore, there exist potive constants C and To such that, for all r > to 

E{t) < C'(J(t) + e _Tr ). 

From the case (i) of Lemma 2.3 one can conclude that e~ rT /J(t) is decreasing as t —^ oo 
provided T > 2dist ( D,B ). Thus we have E{r) = 0(J(t)) as t —»• oo. Applying this together 
with the trivial estimate J(t) = 0(1) as t —> oo to the right-hand side on (30), we finally 
obtain (24). 


2.6 Quantitative state of the surface 

A combination of (25) and the second order term of the asymptotic expansion of J(t) as t —> oo 
yields information about the value of (3 at all the first reflection points, going from the centre to 
dD. 

Theorem 2.3([24]). Let 7 = 0. Assume that dD is C 5 and f3 € C 2 (dD); Aqd(p) is finite 
and satisfies (20). For each q G A qd{p) let e 3 , j = 1,2 be an orthonormal basis of the tangent 
space at q of dD with e\ x e 2 = n q . Choose an open ball U centred at q with radius r q in such 
a way that there exist a h G Cq(R 2 ) with h( 0,0) = 0 and V/i(0,0) = 0 such that U D dD = 
{q + oyei + a 2 e 2 + h(a 1 ,a 2 )n q \ a\ + cj\ + h{ai,a 2 ) 2 < r 2 }. 

If T > 2dist(D, B), then we have 


lim t 5 


A dD{v) ) 

7T77 . , 7T o _ , . 

= od(p) + A 


(31) 


11 



where 


B 9 d(p) = J2 


CdD(q) 


C 9 D{q ) = - 


qeA dD ( P ) y/det(S q (dB dgD{p) (p)) - S q (dD )) 
1 11 — 12 dg D (p)H dD (q) 


+ 


d 9 D(p) 3 8 d aD {pf det{S q {dB dgD ( p) (p)) - S q (dD)) 


and 


^d 9 £) (p) 2 h’ a p a i crr (^)^ l(T s er t (T » (0) 


- Bpg Bq r Bf U ' „ BpgBnf B 


^-ps^qr^-vu, • g- 



+ 16 d dD (p) 2 H °W' (°) S pr fi 9« 


/3(g) 

ddD(p ) 2 


(Bpq) 


\ddD{p) 


h 



Note that we have used the summation convention where repeated indicies are to be summed 
from 1 to 2. The explicit second-order term of an expansion of Laplace type integral J(r) is 
essential and it is an application of an expansion formula of the Laplace type integral. For this 
see [2]. 

As a corollary of Theorem 2.3 we obtain a procedure for calculating the value of ft at a 
known point on A qd{p) provided y = 0. More precisely, we assume that: 

(i) we know in advance a point q £ A 9 d(p)', 

(ii) we know that dD near q is given by making a rotation around the normal at q of a graph 
of a function h defined on the tangent plane at q of dD and that, in an appropriate orthogonal 
coordinates on the tangent plane, say a = ( 07 , < 72 ), the Taylor expansion of the function at a = 0 
has the form h{a\, a?) = X) 2 <|a |<4 h a o a + ■ ■ ■ with known coefficients h a for 2 < |a| < 4. 

Note that, from (i) we know also dg^ip) = |p — q\, n q = (p — q)/\p — q\ and the tangent plane 
(x — q) ■ v q = 0 at q of dD. 

Fix s G]0, dgo(p){- Choose an open ball B' centred at p — sv q and satisfying B' C B flgD ( p \(p). 
Let T > 2dist (D, B'). Generate the wave up by ./ = xb> and observe the wave on B' over time 
interval ]0, T[. Since A go{p ~ sn q ) = {q} and (20) for p replaced with p — sv q is satisfied, one 
gets (31) in which ball B is replaced with B' and p replaced with p — su q . Therefore, we obtain 
Cg D (q) which yields a linear equation with unknown j3(q) and thus solving this, one obtains 
P(q). Note that, in this procedure we do not assume that A goip) is finite. 

Now it is natural to consider the following problem. 

Open problem 2.1. Asuume that, say, 7 is sufficiently smooth on dD. Find a formula for 
calculating 7 at a known q G A 3d(p) from Uf on B'x ]0, T[ generated by / = Xb’i where B' is 
the same as above. 

The point is to find the asymptotic profile of E{r) in (12) as r —> 00 in terms of v on D. For 
one-space dimensional case we have an explicit formula. See [20]. 


2.7 Other wave equations 

In [20] a result analogous to Theorem 2.1 has been established also for the equation 

a(x)d 2 u — A u = 0 inR 3 x ]0, T[ 

provided: 7 > a(x) > a - e - x £ R for a positive constant 7 and a(x) = 1 a.e. x £ R 3 \ D; 
there exists a positive constant C such that a(x) < 1 — C a.e. x £ D or a(x) > 1 + C a.e. 
x £ D. 
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In [16] the equation 


dfu — V • A(x)Vu = 0 inRrx ]0, T[ 

with a 3 x 3 uniformly positive definite real symmetric matrix-valued function coefficient A{x ) 
satisfying A{x) = I 3 a.e. x £ R 3 \Z? has been studied. It is assumed that each component of A(x) 
is essentially bounded and there exists a positive constant C such that (A(x) — Is)£ ■ £ > C |£| 2 
a.e. x £ D and all £ £ R 3 or —(A(x) — /s)£ • £ > (7|£| 2 a.e. x £ D and all £ £ R 3 . Then, it is 
clear that Theorem 1.2 in [16] for this equation yields also a result analogous to Theorem 2.1. 
However, for both equations there is no result corresponding to Theorems 2.2 and 2.3 via the 
Enclosure Method. 


2.8 Interior problem in time domain 

Let D be a bounded domain of R 3 with C 2 -boundary. Given / £ L 2 (D) satisfying supp / C D 
denote by Uf the solution of the following initial boundary value problem for the wave equation: 


d 2 u — Au = 0 in Z?x]0, T[, 

Du 

— —-— j3(x)u = 0 on dDx]0,T[, 

u(x, 0 ) = 0 in D , 

d t u(x, 0) = f(x) in D, 


(32) 


where v denotes the unit outward normal to D on dD, (5 and 7 are the same as those in (1). 
Problem 2.2. Let B be an open ball and satisfy B C D. Generate u = uj by the initial data 
/ = xb and observe u on B over time interval ]0, T[. Extract information about the geometry 
of dD, 7 and f3 from the observed data. 

In [21], we considered the case when 7 = /3 = 0 in (32) and obtained two theorems corre¬ 
sponding to Theorems 2.1 and 2.2. It will be possible to obtain a theorem corresponding to 
Theorem 2.1 for general 7 and /3; theorems corresponding to Theorems 2.2 and 2.3 for 7 = 0 
and general (3. Thus, a real problem to be solved should be the same as Open problem 2.1. 


3 Further applications and problems 

3.1 Bistatic data, spheroid and simultaneous rotation 

Let 0 < T < 00 . Let / £ L 2 (R 3 ) satisfy supp / fl D = 0. Let u = Uf(x,t ) be the solution of the 
initial boundary value problem: 

d 2 u — Au = 0 in (R 3 \ D)x ]0, T[, 

u = 0 on 8 Dx ]0, T[, 

u(x, 0) = 0 in R 3 \ D, 

dtu(x, 0) = f(x) in R 3 \ Zb 

Problem 3.1. Let B and B' be two known open balls centred at p £ R 3 and p r £ R 3 with radii 
7 ] and 7 ', respectively such that B n D = 0 and B D D = 0. Let \B denote the characteristic 
function of B and set / = \B■ Assume that D is unknown. Extract information about the 
location and shape of D from the data Uf(x,t ) given at all x £ B' and t £ ]0, T[. 
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Let xb' denote the characteristic function of B' and set g = xb 1 ■ The results of this 
subsection are concerned with the asymptotic behaviour of the bistatic indicator function Ib,b' 
defined by 

Ib,b’(t)= _(fv g ~ w f g)dx, 

JR 3 \D 

where Wf is the same as (4) and v g is the solution of (2) in which / is replaced with g. Note 
that Ib,b'(b ) can be computed from Wf on B and thus from uj on Bx ]0, T[. 

Define cj)(x;y,y') = \y — x\ + \x — y'\, (x,y,y r ) £ R 3 x R 3 x R 3 . This is the length of the 
broken path connecting y to x and x to y'. We denote the convex hull of the set F C R 3 by [i 7 ]. 

Theorem 3.1([23]). Let [B U B'] n dD = 0 and T satisfy 


T > 


min 

xedD,yGdB,y'edB' 


<!>{x',y,y')- 


(33) 


Then, there exists a tq > 0 such that Ib.b' ( T ) > 0 for all r > tq and the formula 


lim -log Ib,b'(t) 

t - >00 7 - 


min 

x£dD, y£dB, y'£dB 


,(/>{x;y,y'), 


(34) 


is valid. 

It is easy to see that 


J/Y) = m in 4>{x-,p,p') - (y + y')- 

x£oD, y£oB, y'^oB' x£oD 

Thus formula (34) enables us to extract iriin xei 9 D 4>(x;p,p') from Uf(x,t) given at all x £ B' and 
t €]0, T[. 

The quantity min xe g£) f>[x\p,p') coincides with the shortest length of the broken paths con¬ 
necting p to a point q on dD and q to p', that is, the first reflection distance between p and p' by 
D. Thus, Theorem 3.1 yields a mathematical method for extracting the first reflection distance 
from the waveform of the observed wave. 

Given c > \p — p'\ define E c (p,p') = {x £ R 3 | (j)(x;p,p r ) = c}. This is a spheroid with 
focal points p and p'. Given direction uo £ S 2 at p' let C t (u)- 1 p,p') denote the unique point on 
E c {p,p') PI {p' + sui | s > 0}. f(u]p,p,p') has the expression C(w;p,p / ) = p' + s(uj;p,p')uj with a 
unique s(pjyp,p') > 0 and the map S 2 : to i —> ((u]p,p') £ E c (p.p r ) is bijective. 

Define Aqd{p,p') = {q £ dD | min c p(x;p,p’) = cj>(q',p,p')}. One can write Aqb{p,p') = 

xGdD 

dD fi E c (p,p'), where c = mm. x ^Q£, <f{x\p,p'). 

Similary to Section 2.2, using Theorem 3.1, one can make a decision whether given direction 
u> £ S 2 at p' C(u\p,p') which is a point on E c (p,p') belongs to dD or not. It is based on the 
following chracterization of A dD(p,p')- 

Lemma 3.1 (Proposition 5.1 in [23]). Fix s €]0, rf[. We have: 

(i) if Qiy^'iPiP') belongs to dD, then <f(x-,p,p' + su) = c — s; 

(ii) if C(iv,p,p') does not belong to dD, then m.m x& gB 4>(x-,p,p'+ su) > c — s. 

Therefore, we obtain the following characterization of Aqd(p, p')'. 

A 8 d{p,p') = {C(w;p,p') | min f(x:p,p' + suf) = c - s}. 

x £oD 

The procedure for finding C(w;p,p / ) belonging to dD from a single set of the bistatic data is 
the following: 

Fix a large T and s e]0, r/[. 


14 



Step 1. Generate uj by the initial data / = xb and observe Uf on B' over time interval ]0, T[. 
Step 2. Choose an open ball B" C B' centred at p' + su. 

Step 3. Determine vajm x&9 £, (f{x',p,p' + sw) from the restriction of Uf in the first step onto 
B"x ]0, T[ via Theorem 3.1. 

From the computed value min X £.qd 4>(x;p,p'+su) in the third step, one has: if min x ^qd <t>{x',p,p'+ 
slu) = c — s, then ((oj;p,p') belongs to dD ; if not, then £(w;p,p / ) does not belong to dD. 

Therefore, in principle, one can determine all the points in A qd(p,p') from Uf on B'x ]0, T[ 
for / = xb- This is an advantage of the bistatic data not being seen in the monostatic data. 

The next theoretical result is concerned with obtaining information about shape of D. In 
the following theorem, for simplicity of description, we assume that D is convex. In this case 
A8 d(PiP') consists of a single point q = q(p,p'). 

Theorem 3.2([23]). Assume that dD is C 3 . Let c = min xe gD (j)(x;p,p'). Let T satisfy (33). 
Then, we have 


det{S q(py) {E c {p,p')) - S q(py) (dD)) > 0 


and the formida 


lirn T 4 e rmil heaD, !/6 as, !/ 'e9s' 4>(x;y,y') Ibb ,^ 


7 r / diam B\ / diam B'\ 1 

2\2\q-p\) \2\q-p'\) ^ det{S q (E c {p,p>)) - S q {dD)) 


(35) 


(36) 


is valid. 

Remark 3.1. For more general condition on D instead of it’s convexity see [23]. In that case, 
instead of (35) we have to assume that 


det (S q (E c (p,p')) - S q (dD)) > 0 \/q G A 9D (p,p'). 


Note that in that case A qd{p,p') does not necessary consists of a single point and thus its 
finiteness should be assumed. As a result the right-hand side on (36) should be changed. 

As a consequence of Theorem 3.2 we obtain two procedures for extracting the curvatures 
and principle directions. 

First we describe a procedure for extracting the curvatures of dD at q(p,p') provided q = 
q{p,p') is known. Formula (36) gives the information det (S q (E c (p,p > )) — S q (dD)) from Uf on 
B’x ]0, T[. Then by restricting Uf on B”x ]0, T[ we obtain also from (36) det ( S q (E c _ s (p,p ' + 
sw)) — Sq(dD)) for uj = (q — p')/\q — p'\. Using an analogous equation to (22) (Lemma 5.1 
in [23]), we see that these two quantities construct a linear system with the Gauss curvature 
K dD {q) and a modification of the mean curvature of dD at q, that is 


HdD(q;p,p') = H dD (q ) 


Sg{dD)(A q (p) x A q (p')) • (A q (p) x A q (p')) 
2(1 + A q (p) ■ A q (p')) 


where A q (x) = (q — x)/\q — x\. The system corresponds to (23) and always uniquely solvable. 
Thus, by solving the system we obtain those two curvatures. Therefore, one can obtain an 
approximate shape of dD around q(p,p'). 

Second we show that it is posible to obtain the principle curvature directions of dD at q(p,p') 
by making a rotation of B and B' at the same time around the normal at q(p,p'). We denote 
by p(9) and p'{9) the points rotated around the line directed u q at q = q(p,p') counterclockwise 
with rotation angle 0 G [0, 27t[ of p and p'. Thus p(0) =pan p\ 0) = p'. Then, for all 9 £ [0, 27t[ 
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we know that A dD (p(9),p'{9)) = {q{p,p')}, A q (p{6)) ■ A q (p'{9)), \A q {p(9)) x A q (p'(9))\ and 
| p(9) — q\ + \q — p'(9)\ at q = q(p,pl) are invariant with respect to 9. 

Let B{9) = {x € R 3 | \x — p{9)\ < 77 } and B( 0) = B ; B'{9) = {x € R 3 | \x —p'{9)\ < rf} and 
B'{ 0) = B'. We have [ B(9) U B'{9)\ fl D = 0 provided [B U B'] fl D = 0. Let f(9) denote the 
characteristic function of B(9). 

Then, from Uf^ on B'(9)x ]0, T[ we obtain the function of 9: 

H dD (q-,p(9),p'(9 )) = H dD {q) - 1 ~ Aq{p ^' S q (dD)(V{9)) ■ V(9) 

where V(9) denotes the unit vector directed to A q (p(9 )) x A q (p'{9)). 

Now assume that A q (p) x A q {p') / 0. Then, V{9) attains all the tangent vector at q of dD 
and thus from the behaviour of Hgu(q; p(9 ), p'(9)) as a function of 9 one can determine all the 
directions of principle curvatures say, V(9\) and V{92) with some 6 \ and 9 -2. Since we have 

S q (dD)(V(9 1 )) • V{0{) + S q (dD)(V(9 2 )) ■ V(9 2 ) = 2 H dD (q), 

the arithmetic mean of Hgf)(q-,p(9i),p'(9i)) and Hgjg(q;p(9 2 ),p'(9 2 )) coincides with 

H SDiq) . 

Thus we obtain Hgjg(q). Therefore, we can extract S q (dD) from Uf^ over B'(9)x ]0, T[ given 
at all 9 £ [0, 27t[. This is an advantage of the data collection using a simultaneous rotation of 
the emitter and the receiver. 

Open problem 3.1. Extend the results to other boundary conditions, transmission conditions 
(see Section 2.7) or the Maxwell system (see Section 3.2). 


3.2 The Maxwell system 


In this section we briefly comment on a recent application [25] of the Enclosure Method to an 
inverse obstacle problem whose govering equation is given by the Maxwell system in the time 
domain. 

Let 0 < T < 00 . We denote by E and H the elctric and magnetic fields, respectively. 
Assume that E and H are induced only by the current density J at t = 0 and that the obstacle 
is a perfect conductor placed in the whole space R 3 . The governing equations of E and H take 
the form 


dE 

e— - V xH = J in 
dt 

dH 

a— -b V x E = 0 in 

at 


(R 3 \D)x ]0, T[, 
(R 3 \ ~D)x ]0, T[, 


v x E = 0 


on dDx ]0, T[, 


E\ t =o = 0 in R 3 \ D, 


{ H | t=0 = 0 in R 3 \ D, 


where u denotes the unit outward normal to D on dD ; e and p denote the electric permittivity 
and magnetic permeability assumed to be positive constants. 

There are several choices of the current density J as a model of the antenna. In [25] it is 
assumed that J takes the form 


= f(t)xB(x)a, 


(37) 
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where a / 0 is a constant unit vector, \b denote the characteristic function of B and / € 
H 1 ( 0, T) with /(0) = 0; B is an open ball with very small radius and satisfies B n D = 0. 

In [25] the author considered the following problem. 

Problem 3.2. Fix T. Generate E and H by the source J given by (37) and observe E on B 
over time interval ]0, T[. Extract information about the geometry of D from the observed data. 

Two theorems corresponding to Theorems 2.1 and 2.2 have been obtained in [25]. The main 
difference from the scalar case is the existence of directivity of the source at t = 0 and one of 
two theorems catches the effect of the source directivity. 

The boundary condition imposed on the surface of the obstacle is a typical one like the 
Dirichlet boundary condition for the wave equation. As a next step it is natural to ask: how 
about the case when the electromagnetic wave satisfies a more general boundary condition like 
the Leontovich condition on the surface of the obstacle (see , e.g., [1])? 

Open problem 3.2. Consider the Leontovich boundary condition instead of the perfect con¬ 
ductivity condition: 

v x H — X(x) v x {E x u) = 0 on dDx ]0, T[. 

Extract information about the geometry of D and A from the observed data. 
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